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ABSTRACT 

We use a set of Monte Carlo simulations to follow the cascade produced by a primary 
electron of energy E in in the intergalactic medium. We choose E in = 3 — 10 keV 
as expected from the decay of one of the most popular Warm Dark Matter (WDM) 
candidates, sterile neutrinos. Our simulation takes into account processes previously 
neglected such as free-free interactions with ions and recombinations and uses the 
best available cross sections for collisional ionizations and excitations with H and He 
and for electron-electron collisions. We precisely derive the fraction of the primary 
electron energy that heats the gas, ionizes atoms and produces line and continuum 
photons as a function of the ionization fraction. Handy fitting formulae for all the 
above energy depositions are provided. By keeping track of the individual photons we 
can distinguish between photons in the Lya resonance and those with energy E < 10.2 
eV that do not interact further with gas. This separation is important because a Lycv 
background can heat or cool the gas depending on the nature of the photons, and can 
have effects on the 21 cm radiation emitted by neutral H, which will probably become 
detectable at z > 6 in the near future by the next generation radio interferometers. 
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1 INTRODUCTION 

The determination of the gas temperature and ionized frac- 
tion in the intergalactic medium (IGM) and interstellar 
medium (ISM) is fundamental for a number of astrophys- 
ical studies. In particular, it is essential for the investigation 
of the nearly uniform, dark, neutral state of the Universe 
known as Dark Ages, which has become in the last decade 
one of the most studied topics in cosmology. The cosmic 
phase between hydrogen recombination at z ~ 1000 and the 
so-called Epoch of Reionization (EoR) at z ~ 9 can be di- 
rectly detected by the 21 cm hyperfine triplet-singlet level 
transition of the ground state of neutral hydrogen. A new 
generation of low frequency radio interferometers such as 
the Low frequency Array (LOFAR), , the 21 Centimeter Ar- 
ray (21CMA), the Mileura Wide-field Array (MWA) and the 
Square Kilometer Array (SKA), are expected to reach the 
sensitivity required to map the HI distribution at angular 
resolution of the order of a few arcminutes (e.g. Pen, Wu, & 
Peterson 2004; Bowman, Morales, & Hewitt 2005; Kassim et 
al. 2004; Wyithe, Loeb, & Barnes 2005). If a few issues such 
as foreground removal, interferences from human-generated 
signals and ionospheric corrections will be successfully taken 
care of these large radio-arrays will be able to perform a to- 
mography of the Universe before and during the EoR. 

It is therefore crucial to be able to predict the magni- 
tude of the signal. The standard scenario predicts the gas 



temperature, Tk, to decouple from the Cosmic Microwave 
Background (CMB) temperature, Tcmb, at z ~ 300. Tk 
then starts decreasing adiabatically until the first sources 
of light heat again the gas above Tcmb at z ~ 20. In this 
case the 21 cm radiation would be seen in absorption against 
the CMB in the redshift interval 20 < z < 300 as discussed 
in detail by many recent studies (e.g. Loeb & Zaldarriaga 
2004). 

A number of more complex theoretical scenarios have 
been investigated recently. In particular if there was a mech- 
anism to heat the gas at z > 10, thus producing a thermal 
history different from the standard one, a direct imprint 
would be left on the 21 cm radiation from the neutral gas, 
which could then be observable in emission rather than in 
absorption. 

First stars (Ciardi & Salvaterra 2007), intermediate 
mass black holes (Zaroubi et al. 2006; Ripamonti, private 
communication) and decaying or annihilating Warm and 
Cold Dark Matter (CDM) particles (Shchekinov & Vasiliev 
2006; Furlanetto, Oh, & Pierpaoli 2006; Valdes et al. 2007) 
could be capable of leaving a significant trace in the evolu- 
tion of the neutral intergalactic gas. Future 21 cm observa- 
tion will possibly rule out or confirm these theoretical pre- 
dictions, opening a new exciting frontier in cosmology and 
possibly unveiling the nature of dark matter particles. 

All the most popular dark matter particle candidates 
inject energy into the IGM, either via decays or annihila- 
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tions, initiating an energy cascade from energetic primary 
photons or electrons. The energy deposition depends on the 
large number of interactions taking place during the propa- 
gation of the cascade particles through the IGM. It is there- 
fore very important for achieving correct results to follow 
in detail these cascades and to know exactly how much of 
the energy injected ionizes the gas, produces radiation by 
collisional excitations, and heats surrounding medium, re- 
spectively. 

This fundamental problem has received attention in the 
past (e.g. Bergeron & Collin-Souffrin 1973), and the results 
achieved were corrected by an extensive study by Shull and 
van Steemberg (1985), which we will denote hereafter as 
SVS85, which was a development of a previous work by Shull 
(1979), hereafter S79. 

In this work we present the results of a Monte Carlo 
calculation in the spirit of the one from S79 and SVS85; 
in comparison, we use more recent and accurate cross sec- 
tions for collisional ionization and excitations from electron 
impacts, for electron-electron collisions, for free-free inter- 
actions and recombinations. In addition to this we follow in 
detail the radiation produced by the excitations and we are 
able to predict precisely how much of the energy goes into 
photons that do not further interact with the gas (E< 10.2 
eV) and how much of it contributes to the Lya background, 
which directly affects the physics of the 21 cm line radiation 
by the so called Wouthuysen-Field effect (e.g. Wouthuysen 
1952; Field 1959; Hirata 2005) and which can heat or cool 
the gas depending on the nature of the Lya photons (Chen 
& Miralda Escude 2004, Chuzhoy & Shapiro 2007). 



2 METHOD 

2.1 Description of the calculation 

We calculate the effects produced by an X-photon of ~keV 
energy injected into the IGM with T < 10 4 K. At these 
energies the dominant interaction is photoionization (see e.g. 
Zdziarski & Svensson 1989) and the X-photon ionizes an H 
or He atom producing an energetic primary electron. We 
follow the subsequent secondary cascade products. 

For our calculations we choose two specific energies for 
the primary electron, E in = 3 keV and 10 keV. This energy 
range is of great interest because sterile neutrinos, one of 
the most promising WDM candidates, are expected to emit 
line radiation at an energy between 3-25 keV. A large effort 
has been done recently by several authors to constrain the 
mass and lifetime of radiatively decaying sterile neutrinos 
from X-ray observations (Abazajian, Fuller & Tucker 2001; 
Abazajian 2006; Abazajian & Koushiappas 2006; Boyarsky 
et al. 2006b; Watson et al. 2006; Mapelli & Ferrara 2005; 
Boyarsky et al. 2006a). 

Once the primary electron is injected into the IGM 
the code calculates the cross sections relative to a list of 
possible processes: H, He, Hel ionization; H, He excitation; 
collisions with thermal electrons; free-free interactions with 
ionized atoms; recombinations. We then estimate the mean 
free paths and derive the probability for a single electron to 
have any of the aforementioned interactions. We assume, in 
agreement with the 3-yr Wilkinson Microwave Anisotropy 
Probe (WMAP) data analysis, that the helium fraction by 



mass is /jy e = 0.248 (Spergel et al. 2007) and we follow 
the fate of a primary electron and its secondary products 
for different values of the ionized fraction x e of the gas, 
which, as showed by SVS85 and S79, is the free param- 
eter which affects the results. We assume, similarly, that 
x e = n(H+)/n(H) = n(He+)/n(He). 

For each primary energy E in and for each assumed value 
of the gas ionized fraction x e we performed 1000 Monte 
Carlo realizations, a number sufficient to produce consistent 
results not biased by the random nature of the computation. 
We will return on this aspect at the end of this Sec. 

If the primary electron (or an energetic secondary elec- 
tron) collisionally ionizes an H or He atom then the resulting 
two electrons have to be followed separately as they interact 
further with the gas. 

Once the electron energy has degraded below 10.2 eV 
we assume that its entire energy is deposited as heat. This 
is a simplification as the electron would actually thermalize 
with the gas and as a consequence the precise heating should 
be calculated by taking into account the gas temperature. 
If the temperature is of the order of 10 4 K electrons with 
energies lower than 1 eV could even cool the gas. 

Collisional excitations of H ad He produce photons that 
escape freely in the surrounding medium if their energy is 
lower than 10.2 eV or that can further interact with the gas 
if they have higher energy. 

The study from S79 and SVS85 derives the amount of 
energy which is deposited in excitations but does not give 
details about the individual photons. We want to estimate 
instead the amount of energy that goes into Lya photons: 
this radiation can in fact interact further with the gas by the 
Wouthuysen-Field effect and is crucial to correctly estimate 
the 21 cm signal. Furthermore, Lya photons by scattering 
resonantly off neutral hydrogen can cool or heat the gas, 
depending on whether they enter the resonance from its red 
or blue wing respectively, as we will explain more in detail 
later in this Letter. 

An additional novel feature of our model is the inclusion 
of the previously neglected two-photon forbidden transition 
2s — ► Is. We included both the direct collisional excita- 
tion cross section to the 2s level and the probability that a 
collisional excitation to a level n > 3 results in a cascade 
through the 2s level rather than through 2p. As noted re- 
cently (Hirata 2005, Chuzhoy & Shapiro 2007) this effect 
is not negligible in general and constitutes the most proba- 
ble decay channel for H atoms excited to a level n = 3 or 
higher. The result is the emission of two photons below the 
Lya energy that do not further interact with the gas. 

With our code we are therefore able to separate the 
Lya radiation from the less energetic photons produced by 
cascade from the 2s level. All Lya photons that result from 
collisional excitations (or injected Lya photons, see Chen & 
Miralda Escude 2004) have a cooling effect on the gas. 

We also included in our calculation processes that can 
produce continuum photons, such as recombinations and 
Bremsstrahlung free-free interactions of electrons with ion- 
ized atoms. This effects are negligible as we will see and there 
is virtually no production of photons between the Lyman- 
a and the Lyman-/3 resonances. This radiation would have 
redshifted into continuum Lya photons (Chen & Miralda 
Escude 2004) and by entering the Lya resonance from its 
blue wing would have heated the gas instead. 
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Gas 


Excitations 


Ionizations 


Excitations 


(ionized fraction) 


Heating 


(Lyman-a) 


(H, He, Hell) 


(E < 10.2 eV) 



0.00010 
0.00015 
0.00020 
0.00030 
0.00050 
0.00100 
0.00150 
0.00200 
0.00300 
0.00500 
0.01000 
0.01500 
0.02000 
0.03000 
0.05000 
0.10000 
0.15000 
0.20000 
0.30000 
0.50000 
0.70000 
0.80000 
0.90000 
0.95000 
0.99000 



0.136566±0. 
0.143090±0. 
0.148240±0. 
0.155762±0. 
0.166903±0. 
0.183868±0. 
0.196387±0. 
0.204938±0. 
0.219333±0. 
0.239705±0. 
0.271348±0. 
0.295171±0. 
0.314125±0. 
0.343381±0. 
0.388884±0. 
0.470896±0. 
0.535092±0. 
0.590440±0. 
0.676904±0. 
0.807223±0. 
0.900674±0. 
0.937668±0. 
0.970899±0. 
0.985733±0. 
0.997217±0. 



006526 
007234 
006846 
007440 
007935 
008404 
008972 
009070 
009728 
010887 
012790 
013789 
015099 
016290 
019726 
025642 
027713 
029182 
027731 
022754 
014926 
010710 
006796 
004310 
001808 



0.334972±0. 
0.330344±0. 
0.326795±0. 
0.322317±0. 
0.315132±0. 
0.305810±0. 
0.298643±0. 
0.294879±0. 
0.286883±0. 
0.277359±0. 
0.264273±0. 
0.254455±0. 
0.246530±0. 
0.235511±0. 
0.217985±0. 
0.188577±0. 
0.165122±0. 
0.145668±0. 
0.114286±0. 
0.067655±0. 
0.034983±0. 
0.022106±0. 
0.010273±0. 
0.004954±0. 
0.000993±0. 



012904 
012740 
012202 
012472 
012997 
012733 
012996 
012814 
012733 
013153 
013665 
013040 
013916 
013593 
014083 
014297 
014062 
013954 
012867 
009607 
006577 
004776 
002993 
002056 
000824 



0.379666±0. 
0.379530±0 
0.379611±0. 
0.378449±0. 
0.377684±0. 
0.374306±0. 
0.372148±0. 
0.368847±0. 
0.366153±0 
0.359742±0. 
0.347151±0. 
0.337623±0. 
0.330200±0. 
0.316941±0. 
0.296736±0. 
0.257177±0. 
0.226825±0. 
0.199731±0. 
0.158192±0. 
0.095154±0. 
0.048889±0. 
0.030452±0. 
0.014284±0 
0.007120±0. 
0.001352±0 



013697 
012823 
013011 
012658 
013238 
012698 
012737 
012509 
012303 
012514 
012160 
011888 
012187 
012406 
013365 
014095 
016066 
016245 
015455 
013304 
009166 
006889 
004777 
003168 
001382 



0.148796±0. 
0.147036±0 
0.145354±0 
0.143472±0 
0.140281±0. 
0.136016±0. 
0.132822±0 
0.131336±0. 
0.127631±0. 
0.123194±0 
0.117228±0 
0.112751±0 
0.109145±0. 
0.104166±0. 
0.096396±0 
0.083350±0 
0.072962±0 
0.064162±0. 
0.050618±0. 
0.029968±0 
0.015454±0. 
0.009774±0. 
0.004543±0. 
0.002193±0. 
0.000439±0 



006668 
006585 
006477 
006310 
006595 
006497 
006566 
006457 
006581 
006770 
006788 
006453 
006836 
006984 
007000 
006941 
006617 
006753 
006023 
004703 
003125 
002281 
001432 
000976 
000424 



Table 1. Fraction of the energy Ei n of a 10 keV primary electron that is deposited into heat, ionizations, Lyoj line radiation and photons 
with energy E < 10.2 eV after two-photon decay from the hydrogen 2s level. For every value of the ionized fraction 1000 realizations 
were performed. 



Electron-electron collisions between secondary and 
thermal electrons were implemented as in SVS85 and S79 
and we treated similarly the energy distribution of secondary 
electrons following collisional ionization of H and He, so we 
refer the reader to those works for a detailed explanation. 

As we mentioned earlier, we found that 1000 Monte 
Carlo realizations gave consistent results. Going from 50 to 
1000 realizations changed the averaged values by less than 
5% and the respective a by less than 10%. We deemed 1000 
to be a sufficient number of realizations for a stable result. 
The la values of the energy deposition fractions are included 
in Table 1 and they vary from 2.9% to 0.04% of the total 
amount of energy. 



2.2 Cross sections 

For our purposes it is important to use the best available 
cross sections for the several interaction channels considered 
in our calculation. We will give here the references rather 
than entering in the description of physical details that are 
not the scope of this Letter. The cross sections (<Ti) for col- 
lisional ionization of H, He, He+ were taken from Kim & 
Rudd (1994), Shah et al. (1987), Shah et al. (1988). A sim- 
ple functional fit is 



<Ji{x) 



4-Kao 



aln(:r) + b 



K) 



+ c 



Int 



t + 1 



(1) 



where x = E kin /E B is the ratio between the incoming elec- 
tron energy and the binding energy of the atomic electron. 



The collisional excitation cross sections of H and He 
are from Kim & Desclaux (2002) , while for the excitation to 
the 2s level of H we used the work from Bransden & Noble 
(1976). The cross section for Coulomb collisions between 
electrons is from Spitzer & Scott (1969) as in S79, while the 
free- free cross section (cr//) for electrons interacting with 
protons is given by the Bethe-Heitler quantum-mechanical 
Born approximated result (see e.g. Haug 1997), 



daff _ 16aZ 2 r 2 
dk 



In 



Pi+Pf 
Pi -Pf 



(2) 



where m e the electron mass, ro is the classical electron ra- 
dius, a is the fine structure constant, k — E/mc 2 is the 
photon energy in units of m e c 2 , and Pi, Pf are the momenta 
of the incident and scattered electron respectively, again in 
units of m e c 2 . 

As mentioned before we also included the recombination 
cross section oy, to take into account a process which could 
produce continuum photons. Neglecting helium recombina- 
tions, we have that for hydrogen 



*r{v,n) ^3 x 10 10 -^2 cm 2 



(3) 



where gf b is the Gaunt factor of 0(1), v is the emitted radi- 
ation frequency, v e is the electron velocity and n is the level 
at which the electron recombines. 
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3 RESULTS 

Our results are summarized in Table 1 and in Figure 1. In 
the Table we report the fraction of the energy of a 10 keV 
primary electron which, for different values of the gas ionized 
fraction x e , is deposited into heat, Lye* excitations, ioniza- 
tions and photons with E < 10.2 eV. The errors correspond 
the standard deviation from the mean calculated over 1000 
Monte Carlo realizations of the experiment. 

The energy fraction that goes into heating grows rapidly 
as the gas ionization fraction value becomes higher. The 
physical reason of this behavior is that electron-electron in- 
teractions become dominant for high values of x e . 

We performed our calculations for 3 keV and 10 keV 
to account for the theoretically and observationally inferred 
range of the sterile neutrino mass and found that the energy 
fraction distributed among the different processes remains 
similar in both cases (within 2%). This is in agreement with 
the results from S79 and SVS85 that find that for primary 
electrons with energies higher than 100 eV the fractional 
energy depositions converge rapidly to a common behavior. 

We gave a convenient functional form to our results by 
fitting the data in Table 1 with an accuracy greater than 
3.5%. 

• Fraction fh of the primary energy deposited into heat: 
f h = 1.0 - 0.8751 ( 1.0 -a; ' 4052 ) (4) 

• Fraction f a of the primary energy converted into Lya 
radiation: 

f a = 0.3484 ( 1.0 - a; - 3065 ) - 9533 (5) 

• Fraction /» of the primary energy deposited into ioniza- 
tions: 

fi = 0.3846 ( 1.0 _ ^420)1.1952 (6) 



• Fraction f c of the primary energy deposited into con- 
tinuum radiation: 

f c = 0.1537 (1.0 -a: ' 3224 ) (7) 

Our results differ sensibly from those of the past stud- 
ies, with differences as high as 30% for some values of x e . 
The reason for this is a combination of the better and more 
modern cross sections which we used and of the additional 
processes included in the calculation such as the transitions 
through the 2s level of H. We compare in detail our results 
with those obtained by SVS85 in Fig 2. While for small 
values of x e the curves are similar it is evident that the dif- 
ferences become substantial for x e £ > 0.1 



4 DISCUSSION AND CONCLUSIONS 

We have presented an updated calculation of the energy cas- 
cade arising from a primary electron with energy in the range 
(3 — 10 keV) predicted for one of the most popular dark 
matter particle candidates, i.e. sterile neutrinos. We have 
computed the fractional energy deposition into ionizations, 
excitations and heating to a new level of detail and followed 
the fate of individual photons to be able to distinguish be- 
tween Lya injected radiation and continuum photons with 
energies under 10.2 eV. 

As mentioned previously in this work, Lya radiation 
affects the gas in several ways. It has a thermal effect on the 
matter, with a heating or cooling effect depending on the 
nature of the photons. Line or injected Lya photons cool 
the gas, while photons between the Lya and Ly/3 resonances 
that redshift to 10.2 eV (continuum Lya photons) act as a 
heat source for the gas. 

A Lya background is also responsible for the 
Wouthuysen-Field process, which directly affects the physics 
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Figure 2. Comparison of the present results (solid lines) with SVS85 (dotted) in terms of the fractional energy of the primary electron 
going into heat, fh, ionizations, /j, and excitations, f e . In the plot we indicate with * the fractional energy depositions as fitted in SVS85. 



of the 21 cm line radiation from neutral H and can make it 
visible in emission or absorption against the CMB. This as- 
pect is important because 21 cm observations of the high 
redshift Universe will be performed in the near future by 
next generation low frequency radio interferometers such as 
LOFAR. 

We included in our calculations the only mechanisms 
that could produce continuum Lya radiation, recombina- 
tions and free-free interactions with ions. Both processes are 
more probable as the ionized fraction increases, but at the 
same time the cross section for electron-electron collisions 
becomes dominant, so we found that both these channels 
are practically negligible and that electrons injected in a 
highly ionized gas are thermalized before they can produce 
continuum photons. 

We expect that the same calculations performed for rel- 
ativistic electrons could produce interesting results in this 
sense, also taking into account that processes such as inverse- 
Compton on CMB photons would become important and 
produce continuum radiation. This could be a useful exten- 
sion to this work and could be applied to study the effects of 
Light Dark Matter decays/annihilations in the energy range 
around ~ 10 MeV. 
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